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On the massless "just-so" solution to the solar neutrino problem 
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We study the effect of the non-resonant, vacuum oscillation-like neutrino flavor conversion in- 
duced by non-standard flavor changing and non-universal flavor diagonal neutrino interactions with 
electrons in the sun. We have found an acceptable fit for the combined analysis for the solar ex- 
periments total rates, the Super-Kamiokande (SK) energy spectrum and zenith angle dependence. 
Phenomenological constraints on non-standard flavor changing and non-universal flavor diagonal 
neutrino interactions are considered. 
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Non-standard neutrino interactions with matter can 
generate neutrino flavor oscillations. This phenomenon 
was suggested by Wolfenstein in his seminal 1978 pa- 
per Applications of this idea to the solar neutrino 
problem were first suggested in 1991 BM when it was 
observed that resonantly enhanced neutrino oscillations 
induced by non-standard neutrino flavor changing (FC) 
as well as non-universal flavor diagonal (FD) neutrino 
interactions can explain the solar neutrino experimen- 
tal data B which clearly indicates a solar neutrino flux 
smaller than what is predicted by the standard solar 
models M. 

Interestingly enough, such oscillations can be reso- 
nantly enhanced even if neutrinos are massless and no 
vacuum mixing angle exists 0, as a result of an inter- 
play between the standard electroweak neutrino charged 
currents and non-universal non-standard flavor diagonal 
neutrino interactions with matter. In fact, in this mech- 
anism, resonance plays a crucial role in order to provide 
a viable solution to the solar neutrino problem |i] || . 

It should be emphasized that if such nonstandard neu- 
trino FC and FD interactions exist only with electrons, 
no resonant conversion can happen because the mixing 
angle in matter is constant, as we will see later, contrary 
to the case of usual MSW effect ||, or the case with 
d,w-quarks FC and FD interactions. From this point of 
view, the oscillation induced by non-standard neutrino 
interactions with electrons alone is similar to the vac- 
uum oscillation mechanism despite the difference that it 
occurs only in matter, inside the sun. 

This non-resonant neutrino conversion can be useful 
to explain the solar observations. The first discussion on 
this possibility appeared in Ref. where non-resonant 
neutrino oscillation induced by FC and FD interactions 
only with electron in the solar matter was mentioned as 
a possible solution to the solar neutrino problem. 

Nevertheless, so far, no quantitative analysis of this 
scenario was presented. 

In this brief report we investigate this possibility by 
performing a detailed fit to the most recent solar neu- 



trino data. We conclude that non-resonant neutrino os- 
cillations induced by non-standard neutrino interactions 
can only provide a rather poor fit to the total rates ob- 
served by all the solar neutrino experiments coming from 
Homestake, Gallex/GNO, Sage and SK [Q whereas when 
we include also the full SK recoil electron spectrum and 
the zenith angle dependence the fit become an accept- 
able one. We find also that this fit requires the new 
non-standard neutrino interactions parameters to be very 
large. 

Here we assume that neutrinos have non-standard FC 
as well as FD interactions only with electrons which 
could be realized in some models such as minimal super- 
symmetric standard model without R-parity JlO[ or 
SU(3) C ® SU(3) L ®U(1) N (331) models @. Under this 
assumption, the evolution equation for massless neutri- 
nos in matter can be expressed as j2|, 
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where A e (r) and Ai(r) (I = (m,t) are, respectively, the 
probability amplitudes to detect a v e and vi at position 
r and 
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describe, respectively, the relative strength of the FC 
and the FD (but non- universal) interactions where G VaVfj 
{a, (3 = e,/i, r) denotes the effective coupling of the re- 
spective interaction. 

In this mechanism the mixing angle in matter 9 m does 
not depend on the electron density and is simply given 
by, 



sin 20 m = 



4r 2 



(l-e') 2 +4e 2 ' 



(3) 



We see that no MSW like resonance can occur be- 
cause the mixing angle in matter is constant and does 
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FIG. 1. Definitions of the variables r and (f>. (Size of the 
neutrino production region was enlarged just for the purpose 
of illustration.) 

not change along the neutrino trajectory (however see 
Ref. @). 

Let us introduce the two variables r and 4> which define 
the production point of neutrinos in the sun as shown in 
Fig. 0. Then, for given values of (e, e') and a given pro- 
duction point in the sun defined by r and </> the survival 
probability of electron neutrinos at the solar surface can 
be written as |6|, 



P{v e — * i/ e ; r, 4>) = 1 — sin 2 29 m sin' 



tf(r, 4>) 



(4) 



where 



f(r, 4>) = v/4e 2 + (1 - e') 2 V2G F / N e (r, cf>, x)dx, 



(5) 

where the N e (r, (f>, x) is the electron density profiles along 
the neutrino trajectory which starts at the creation point 
(r, 4>) corresponding to x — and ends at the solar surface 
corresponding to x — x max . Note that there is no energy 
dependence in the probability. 

From Eq. ([|) we can estimate the oscillation length as 
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2.4 x 10 2 
y/(l - e') 2 + 4^ 
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km, 



(6) 



where we take N e = N e (R ~ O.li?©) ~ 65 mol/cc as 
a reference value. From Eq. (|^) we see that if either 
|1 — e'| or |e| is of the order of 0.01, the oscillation length 
is typically less than a few percents of the solar radius 
in the neutrino production region. This implies that for 
such values of e and e' there are many oscillations before 
neutrinos reach the solar surface and the final survival 
probability which is averaged over the neutrino produc- 
tion point will be, 
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FIG. 2. Allowed parameter region. Region allowed by the 
total rates. Best fit is obtained when lei = 2.4 x 10 -4 and 



\<l' - 1| = 2.7 x 10" 4 with xl 



7.5 for 2 d.o.f. 



for any values of r and <p and therefore, for any sources 
of neutrinos ||. Therefore, such rapid oscillation cannot 
fit the solar neutrino data well. 

As pointed out in Ref. H, an interesting possibility 
remains if both |1 — e'\ and |e| is smaller than ~ 0.01. For 
such e' and e, if ^ ~ (2n + 1)tt with small n, neutrinos 
produced as u e can be almost v x (x = /i, r) at the solar 
surface. On the other hand, if 'F ~ 2ri7r v e remains as v e 
at the solar surface. 

Since neutrinos from different nuclear reaction origins 
have different production distributions, there is a possi- 
bility that properly adjusting the parameter e and e' neu- 
trinos from different reaction origins could be oscillated 
into another flavor in such a way that the solar neutrino 
data can be accounted for. In principle, this could hap- 
pen if neutrinos oscillate only once or a few times before 
they reach the solar surface, similar to what happen to 
the case of the vacuum long- wavelength oscillation solu- 
tion to the solar neutrino problem fll3|| . 

In order to settle this issue we have performed a de- 
tailed x 2 analysis using the latest standard solar model 
by Bahcall et al. (| (BP98 SSM) as well the latest re- 
sults of the current solar neutrino experiments coming 
from Homestake, Gallex/GNO, SAGE and SK [|. 

The fitting procedure is as follows. We first compute 
the survival probability of v e at the solar surface by the 
formula in Eq.(Q) for various different production points 
defined by r and <f> as in Fig. Then we compute av- 
eraged survival probabilities for neutrinos from different 
sources, i.e., pp, pep, 7 Be, 8 B, 13 N and 15 (we neglect 
other minor contributions from 17 F and hep for simplic- 
ity) taking into account the neutrino production point 
distribution from BP98 SSM. After we compute the ex- 
pected solar neutrino signal for each experiment, using 
the survival probability obtained above, we perform a x 2 
analysis following the prescription given in Ref. flljj ]. 

We show in Fig. Hthe allowed parameter region deter- 
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mined by our % 2 analysis. We have used only the total 
observed rates of solar neutrinos by four experiments. 
The best fit is obtained at (|e|,|e'-l|) = (2.4, 2.9) x 
with Xmin — 7-5 for 4 (data) -2 (free parameters) = 2 
d.o.f. which corresponds to 2.4 % C.L. indicating a poor 
fit. This is because the integrations of the survival prob- 
ability over the variables r and <fi tend to kill the just-so 
suppressions of the neutrino fluxes and the final averaged 
probabilities from different sources end up with rather 
similar values to each other. 

However, the situation is still better than energy inde- 
pendent suppression since the following relation for the 
final averaged survival probability, 

(P( 8 B)) < (PfBe)) < (P(pp)) , (8) 

holds in this mechanism. In fact at the best fit point, 
we have obtained (P( 8 B)) - 0.42, (P( 7 Be)) ~ 0.46 and 
(P(pp)) ~ 0.57. 

We also performed a x 2 analysis allowing 8 B flux to 
vary freely but we do not find any significant improve- 
ment of the fit. 

Let us to include also the spectrum and zenith angle 
dependence in our fit. First, let us note that this mecha- 
nism does not distort the SK energy spectrum since the 
conversion probability is completely energy independent. 
Therefore, contribution in x 2 from the spectrum is just 
constant (does not depend on e' and e) and the allowed 
parameter region shown in Fig. 2 is not affected by the 
spectrum. Second, for the range of parameters we are 
considering in this work, there would be no significant 
effect of the earth matter because the oscillation length 
in the earth is much larger than the earth radius. So 
again, contribution in \ 2 from the zenith angle depen- 
dence is also just constant and the allowed parameter 
region shown in Fig. 2 is not affected by the zenith angle 
dependence. 

The total combined Xmin can be computed by simply 
adding the two constant contributions from spectrum and 
zenith angle without affecting the allowed parameter re- 
gion presented in Fig. 2. In this case, despite the fact 
that the fit is poor only with the total rates, when we 
combined SK spectrum as well as zenith angle depen- 
dence, we have obtained Xmin = ^ or ^ d-o.f. which 
corresponds to 34.4 % C.L.. 

Here let us consider, as an interesting exercise, that 
the case when the systematic error of the Homestake is 
assumed to be 3 times larger than it has been reported. 
In Fig. 3 we present the region allowed by the rates under 
this assumption. We have obtained Xmin = 3-3, for rates 
only which indicate the significant improved over the case 
presented in Fig. 2. We notice that this kind of exer- 
cise could be worthwhile to consider to take into account 
the possibility of some unknown systematic effect of the 
Homestake experiment as it has not been calibrated with 
a radioactive source. 

The main problem appearing in the solution to the 
solar neutrino problem based on FC and non-universal 
FD neutrino interaction with electrons is related with 
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FIG. 3. Same as in Fig. 2 but with the systematic 
error of the Homestake experiment was assumed to be 3 
times larger. Best fit is obtained when |e| — 2.3 x 10 -4 and 
|e' - 1| = 3.1 x 1CT 4 with Xmin = 3.3 for 2 d.o.f. 



the magnitude of the phenomenologically required non 
standard parameters. Our statistical analysis shows that 
although the FC parameter e does not need to be very 
high, (e e i f» 10~ 3 ), the non-universal FD parameter e' el 
is found to be of order of 1. 

The value for the FC parameter e is compatible with 
the available phenomenological tests to the flavor con- 
servation law. In fact, the most stringent constraints on 
this parameter are due to the upper bounds on fi~~ — > 
e~ e + e~ and r~ — * e~ e + e~ fll5|] : 

BR(/x- -► e- e+ e~) < 1.0 x 10~ 12 , 

BR(t" -> e^e+O < 2.9 x 10~ 6 , (9) 

at 90% C.L.. Normalizing the above bounds to the mea- 
sured rates of the related lepton flavor conserving decays, 
BR(//- ->e-^)» 100 % and BR(r~ -> e~ v e v T ) = 
0.1781 H|, we obtain § 

e ep = G^IG F < 1.0 x 10- 6 , 

e er = G eT /G F < 4.2 x 10~ 3 . (10) 

Note, furthermore, that these bounds on e can be also 
relaxed by a factor of 5-6 due to the breaking of the 
SU(2)l symmetry ||. Therefore, assuming that the neu- 
trino transitions involve the first and the third families, 
the required value of e is compatible with the phenomeno- 
logical limits. 

The challenge to this solution is related with the re- 
quired value of the parameter e' el since universality ex- 
perimental tests in the leptonic sector are very much 
stringent. In Ref. [|l6| the constraints involving the sec- 
ond and third lepton families, i.e., interactions involving 
transitions of the type «-» v T are obtained. It was 
found that JT|], 

e^ T < 3.8 x 10~ 3 . (11) 
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Note, however, that the parameter relevant for our 
present analysis of the solution to the solar neutrino prob- 
lem involves necessarily the first neutrino family (v e )- 
Such constraint can be obtained following the same steps 
of Ref. |(| . No direct limit can be obtained to e' er . Nev- 
ertheless, since e' eT = e' — e' , limits on this parame- 
ter are found considering the experimental constraints of 
Eq. ( pi] ) and limits on e' . 

Non-zero values for e ee (e w ) gives a additional contri- 
bution to v e e — > v e e {y^e — > v^e) cross section and can 



put constraints on e ee and e,, 



We use the more re- 



cent data about the v e e — » v e e total cross section [ |L7| . 
This cross section is a function of e e ; and v JGf (the 
axial part of the effective coupling of the respective in- 
teraction). We obtained 



2.56 < G VeV JG F < 0.63, 



(12) 



at 90% C.L. for arbitrary G^ v /Gf- Taking the value 
quoted by Ref. ||: -0.18 <°G ViiV JG F < 0.14. We 



obtained that e' is bounded to 



0.81 < e' efi < 2.70, 



(13) 



where the limit is at 90% C.L. .Using the constraint from 
Eq. (Ill]) and from Eq. (pi), we get finally, 



1.81 < e' - 1< 1.70, 



(14) 



at 90% C.L. . From this constraint, we conclude that 
is possible to satisfy the experimental constraints of FD 
couplings and at same time to be compatible with the 
allowed region of the solar neutrino analysis. A addi- 
tional caution is necessary because the same FD cou- 
plings that induce neutrino oscillations can also change 
the detection cross section, (<j(y e e — > v e e)) used for SK. 
We check that the absolute values of elastic cross section 
inside the range given in Eq.( |i~2| ) are compatible with 
the assumed theoretical errors used in the solar neutrino 
analysis. Also the shape of recoil electron in SK is not 
changed significantly due the FD couplings. 

Concluding, we showed here for the first time a quanti- 
tative analysis of non-standard flavor changing and non- 
universal flavor diagonal neutrino interactions with elec- 
trons as a possible candidate to solve the solar neu- 
trino problem. If the parameters, \e' el — 1| <C 1 and 
e e i 10~ 4 — 10~ 3 then we can get an acceptable fit 
for the combined analysis of the solar experiments to- 
tal rates, the SK energy spectrum and the SK zenith 
angle dependence. We conclude that the constraints on 
the violation of universality allow us a small e eT and a 
large value for e' eT and is compatible with the preferred 
values of our solar neutrino analysis. In this solution, no 
spectrum distortion, no zenith angle dependence and no 
seasonal effects are expected. Also only negative results 
are expected in long-baseline experiments due the very 
large oscillation length. 
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